Unbound or Distant Planetary Mass Population 
Detected by Gravitational Microlensing 

T. Sumi 1 ' 2 - 3 , K. Kamiya 1 ' 3 , A. Udalski 4 ' 5 , D.P. Bennett 1 ' 6 , LA. Bond 1 ' 7 , F. Abe 1 - 3 , C.S. 
Botzler 1 ' 8 , A. Fukui 1 ' 3 , K. Furusawa 1 ' 3 , J.B. Hearnshaw 1 ' 9 , Y. Itow 1,3 , P.M. Kilmartin 1 ' 10 , 
A. Korpela 1 ' 11 , W. Lin 1 - 7 , C.H. Ling 1 ' 7 , K. Masuda 1 ' 3 , Y. Matsubara 1 - 3 , N. Miyake 1 ' 3 , M. 
Motomura 1,3 , Y. Muraki 1 ' 12 , M. Nagaya 1,3 , S. Nakamura 1,3 , K. Ohnishi 1 ' 13 , T. Okumura 1,3 , 
Y.C. Perrott 1 ' 14 , N. Rattenbury 1 ' 8 , To. Saito 1 ' 15 , T. Sako 1 ' 3 , D.J. Sullivan 1 ' 11 , W.L. Sweatman 1 - 7 , 
P.J. Tristram 1 ' 8 , PCM. Yock 1 ' 8 , M.K. Szymanski 4 - 5 , M. Kubiak 4 ' 5 , G. Pietrzyflski 4 ' 5 - 16 , R. 
Poleski 4,5 , I. Soszynski 4 ' 5 , L. Wyrzykowski 4 ' 17 , K. Ulaczyk 4,5 

1 Microlensing Observations in Astrophysics (MO A) Collaboration. 

2 Department of Earth and Space Science, Osaka University, Osaka 560-0043, Japan. 

3 Solar- Terrestrial Environment Laboratory, Nagoya University, Nagoya 464-8601, Japan. 

4 Optical Gravitational Lens Experiment (OGLE) Collaboration. 

5 Warsaw University Observatory, Al. Ujazdowskie 4, 00-478 Warszawa, Poland. 

6 Department of Physics, University of Notre Dame, Notre Dame, IN 46556, USA. 

7 Institute of Information and Mathematical Sciences, Massey University, Private Bag 102- 
904, North Shore Mail Centre, Auckland 0745, New Zealand. 

8 Department of Physics, University of Auckland, Private Bag 92019, Auckland 1142, New 
Zealand. 

9 Dept. of Physics and Astronomy, University of Canterbury, Christchurch 8140, NewZealand. 

10 Mt. John University Observatory, University of Canterbury, P. 0. Box 56, Lake Tekapo 
8770, New Zealand. 

11 School of Chemical and Physical Sciences, Victoria University, Wellington 6I4O, New 
Zealand. 

12 Department of Physics, Konan University, Nishiokamoto 8-9-1, Kobe 658-8501, Japan. 

13 Nagano National College of Technology, Nagano 381-8550, Japan. 

14 Cavendish Laboratory, Cambridge University, J. J. Thomson Avenue, CB3 OHE Cam- 
bridge, UK. 

15 Tokyo Metropolitan College of Industrial Technology, Tokyo 116-8523, Japan. 

16 Universidad de Concepcion, Departamento de Astronomia, Casilla 160-C, Concepcion, 
Chile. 

17 Institute of Astronomy, Cambridge University, Madingley Rd., CB3 OH A Cambridge, UK. 

Since 1995, more than 500 exoplanets have been detected using different 
techniques 1,2 , of which 11 were detected with gravitational microlensing—^. Most 
of these are gravitationally bound to their host stars. There is some evidence 
of free-floating planetary mass objects in young star-forming regions^&££, but 
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these objects are limited to massive objects of 3 to 15 Jupiter masses with large 
uncertainties in photometric mass estimates and their abundance. Here, we re- 
port the discovery of a population of unbound or distant Jupiter-mass objects, 
which are almost twice (1.8^o 8 ) as common as main-sequence stars, based on 
two years of gravitational microlensing survey observations toward the Galactic 
Bulge. These planetary-mass objects have no host stars that can be detected 
within about ten astronomical units by gravitational microlensing. However a 
comparison with constraints from direct imaging 9 suggests that most of these 
planetary-mass objects are not bound to any host star. An abrupt change in the 
mass function at about a Jupiter mass favours the idea that their formation pro- 
cess is different from that of stars and brown dwarfs. They may have formed in 
proto-planetary disks and subsequently scattered into unbound or very distant 
orbits. 

In a gravitational microlensing event, a foreground lens object is detected as a result 
of the characteristic magnification of a background source star as it passes behind the grav- 
itational field of the lens^. The lens object is detected by means of its mass and not its 
luminosity. The duration of the magnification is parameterized by the Einstein radius cross- 
ing time, £e ~ \J M/Mj days, where Mj = 9.5 x 1O~ 4 M is Jupiters mass. Thus, microlensing 
can detect faint planetary mass objects — which are either unbound to any host star— ^ or 
are in very wide orbits^ — as short-timescale events with £e < 2 days. Although £e also 
depends on the distance and transverse velocity of the lens (see Supplementary Information), 
the observed distribution can be used as a statistical probe of the mass function of the 
lens objects because the spatial and velocity distributions in the Galactic disk and bulge are 
reasonably well known. 

The Microlensing Observations in Astrophysics (MOA)— and Optical Gravitational 
Lensing Experiment (OGLE)— groups both conduct microlensing surveys toward the Galac- 
tic bulge. The second phase of MOA, called MOA-II, carries out a high cadence photometric 
survey of 50 million stars in bulge fields with a cadence of 10-50 minutes. This strategy 
enables MOA to detect very short events with < 2 days, which were quite rare in previous 
microlensing surveys that had lower cadences^^ 1 ^ 1 ^. 

In this analysis of the 2006-2007 MOA-II data set, light curves of genuine microlensing 
events were distinguished from intrinsic variables and artefacts by several empirical criteria, 
which have been developed in previous microlensing surveys^^^^'i&i 9 -. The light curves 
must have a single brightening episode of more than three consecutive significant data points 
with a constant baseline, and should be well fitted by a theoretical microlensing model 10 with 
a well constrained (see the Supplementary Information). 

Although there are a thousand microlensing events in this sample, only 474 well charac- 
terized events have passed our strict selection criteria. Ten of these events have < 2 days 
(see Fig. [1] and Table [I]) thus indicating planetary-mass lenses. We have confirmed that 
this event sample has no significant contamination by possible background effects includ- 
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ing: (1) cosmic-ray hits, (2) fast-moving objects, (3) cataclysmic variables (4) background 
supernovae, (5) binary microlensing events, and (6) microlensing by high-velocity stars and 
Galactic halo stellar remnants. For example, effect (1) is excluded because cosmic-rays never 
hit the same place in four consecutive images, microlensing model fits for effects (2) to (5) 
produce a high \ 2 an d unphysical values of parameters, and effect (6) is excluded by proper- 
motion and radial- velocity observations (see Supplementary Information). After the MOA 
event selection was complete, the MOA group requested additional independent light-curve 
data of these short events from the OGLE group. Seven of the ten events with t E < 2 days 
were also observed by OGLE-III—, and none of them have any other brightening in the eight- 
year OGLE-III light curves. For six of these seven events, there are OGLE data obtained 
during the lensing event that confirm the predictions of the MOA microlensing models. Thus, 
the OGLE data confirms the microlensing interpretation of these short events. 

The detection efficiencies for this analysis were estimated with a Monte Carlo simula- 
tion—^. We simulated 20 million artificial events to evaluate the detection efficiency as a 
function of t E , yielding e(t E ) ~ 1%, 3%, 5%, 10%, 15% and 10% at t E = 0.3, 1, 2, 10, 30 and 
100 days, respectively. The details of the efficiency calculations and consistency tests of the 
selected event distribution are discussed in the Supplementary Information. 

The observed t E distribution is compared to two mass function models in Fig. [2] A 
model t E distribution, @(t E ), can be calculated for an assumed mass function with a standard 
Galactic mass density and velocity model— ^>^. We consider two mass functions. The 
first is a broken power-law— ^ dN/dM = M~ a , with the power index of a\ = 2.0 for 
0.7 < M/M Q < 1, a 2 = 1.3 for 0.08 < M/M & < 0.7 and a 3 as a fitting parameter for 
the brown dwarf regime 0.01 < M/M & < 0.08. The second is a log-normal function 22 
diV/dlogM = exp[(logM — log M c ) 2 /(2cr 2 )] with a mean mass M c and a width in logM of 
<7 C , for < M/Mq < 1.0. For both mass functions, we assume that stars that were initially 
above 1 M & have evolved into stellar remnants — white dwarfs, neutron stars or black holes, 
depending on their initial masses— (see Fig. [2] and Supplementary Table 3). 

The mass functions were constrained by a likelihood analysis, with the likelihood func- 
tion given by the product of the model probability ${t E ) of finding N ohs = 474 events with 
each of the observed t E ^, that is; L = Il^° bs ^>{t E ^)e(t E) i). 

We have evaluated the likelihood distributions for these mass functions both with and 
without the tE < 2 events, but the inclusion of the events with t E < 2 days makes little 
difference. The results are shown in Supplementary Table 3 and Supplementary Figs 6 
and 7. Fig. [2] indicates that both models match the data well for t E > 2 days, but at 
t E < 2 days, the ten observed events are well above the model predictions. The power-law 
and log-normal models predict 1.5 and 2.5 events with t E < 2 days, respectively, and the 
corresponding Poisson probabilities for the ten observed events are 4 x 10~ 6 and 3 x 10~ 4 . 
Thus, we feel confident in adding a new planetary-mass population. 

For simplicity, we chose a planetary-mass function model with a ^-function in mass 
Mpl and a fraction of all objects in the planetary-mass populations $pl- The values 
of (M PL /M ,$ PL ) derived from the likelihood analysis are (l.ltoil x 10~ 3 , 0.49±g;^|) and 
(0. 8315x51 x 10~ 3 , 0.46lo;{5) for the power-law and log-normal models, respectively. The 
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contours are shown in Fig. [3j Both models for <?(£e) provide good fits to the entire £e 
distribution, as shown in Fig. [2j The power-law and log-normal models imply 1.9i ;l an d 
1.8l ;g times as many unbound or distant Jupiter-mass objects as main-sequence stars in 
the mass range 0.08 < M/M < 1.0, respectively. These planetary-mass objects are at 
least 1.5 times as frequent as planets with host stars (see Supplementary Information). We 
tested a third mass function that has fewer massive stars and brown dwarfs^, and found 
that the resultant planetary-mass function parameters are consistent with above values (see 
the Supplementary Information). 

The lenses for these short events could be either free-floating planets or planets with 
wide separations of more than about ten astronomical units (AU) from their host stars, 
for which we cannot detect the host star in the light curves^. However, direct imaging, 
with adaptive optics, of planets orbiting young stars places upper limits on planets at wide 
separations. The Gemini Planet Imager has set upper limits^ on the number of stars with 
Jupiter-mass planets at semi-major axes of 10-500 AU. From these results, we estimate that 
<0.4 of the 1.8 planetary- mass objects per star are likely to be bound to stars at orbital 
separations of <500 AU (see Supplementary Information Section 8). Hence, more than 75% 
of these planetary mass objects are probably unbound to stars if their typical mass is a 
Jupiter-mass or more. 

Since the 5-function planetary models are not likely to be realistic, we also tested a 
fourth mass function that is identical to the first, broken-power-law model except for having 
a power-law form in the planetary mass regime, M < 0.01M Q . This yields a planetary-mass 
index of ctpL = 1.3io;4, which is much steeper than the brown dwarf slope of a 3 = 0. 49^27, 
indicating that they are distinct populations (see Supplementary Information). 

Planet-formation theories predict that dynamical instabilities in planetary systems with 
multiple giant planets could scatter many of these planets into unbound orbits^, as well as 
some into large separations^. Recent observations also indicate that planet-planet scat- 
tering plays an important part in moving giant planets into short-period orbits^22. The 
planetary-mass population that we have identified here may have formed in protoplanetary 
disks at much smaller separations and then been scattered into unbound or very distant 
orbits. 



Supplementary Information accompanies the paper on www.nature.com/nature. 
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Figure 1. — Light curves of event MOA-ip-3 and MOA-ip-10. These have the highest signal 
to noise ratio amongst the ten microlensing events with £e < 2 days (see Supplementary 
Fig. 1 for the others). MOA data are in black and OGLE data are in red with error bars 
indicating the s.e.m. The green lines represent the best-fit microlensing model light curves. 
For each event, the upper panel shows the full two-year light curve, the middle panel is a 
close-up of the light curve peak, and the bottom panel shows the residuals from the best-fit 
model in units of magnification, AA. uo indicates the source-lens impact parameter in unit 
of the Einstein radius. The second phase of MOA, MOA-II, carried out a very high cadence 
photometric survey of 50 million stars in 22 bulge fields (of 2.2 deg 2 each) with a 1.8m 
telescope at Mt. John Observatory in New Zealand. MOA detects 500-600 microlensing 
events with 8 months observation every year. In 2006-2007, MOA observed two central 
bulge fields every 10 minutes, and other bulge fields with a 50 minutes cadence, which result 
~8250 and 1660-2980 images, respectively. This strategy enabled MOA to detect very short 
events with £g < 2 days. Since 2002, the OGLE-III survey has monitored the bulge with 
the 1.3-m Warsaw telescope at Las Campanas Observatory, Chile, with a smaller field-of- 
view but better astronomical seeing than MOA. The OGLE-III observing cadence was 1-2 
observations per night, but the OGLE photometry is usually more precise and fills gaps in 
the MOA light curves due to the difference in longitude. 
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Figure 2. — Observed and theoretical distributions of the event timescale, £e- The black 
histogram represents the number N of observed 474 microlensing events in each bin with 
error bars indicating the s.e.m. The red and blue lines indicate the best-fit models with (1) 
the power-law and (2) log-normal mass functions, respectively. For both mass functions, 
we assume that stars that were initially above 1 M & have evolved into stellar remnants - 
white dwarfs, neutron stars or black holes depending on the initial mass. The number of 
remnants is determined by extending the upper main sequence power-law ot\ = 2.0 to 100 
Mq, and the final remnant mass distributions are given by Gaussians^ (see Supplementary 
Table 3) . Each model is multiplied by the detection efficiencies. In each model, dashed lines 
indicate models for stellar, stellar remnant and brown dwarf populations, and the dotted lines 
represent the planetary-mass population (PL). Solid lines are the sums of these populations, 
and both models fit the data well. 
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Figure 3. — Likelihood contours for the planetary-mass function parameters. $ P l indicates 
the fraction of all objects in the planetary-mass population, not including the brown dwarfs 
that have planetary mass in the tail of the log-normal mass function. M P l represents their 
mass. The two sets of contours indicate the 68% and 95% confidence levels. The red and blue 
curves indicate the power-law and log- normal mass functions, respectively, and " +" indicates 
the maximum likelihood points. The top axis scale is Jupiter masses, MJ = 9.5 x 1O _4 M , 
while bottom axis scale is in Solar masses, M & . For the power-law model, the likelihoods 
are evaluated in the (a 3 , M PL , $pl) space and projected into the (M PL , $ PL ) plane. The 
M c = 0.12 and a c = 0.76 parameters are fixed for the log-normal model. The median and 68% 
confidence intervals of (M PL /M Q , <Z> PL ) are (l.lloJ x 10~ 3 , 0.49±g;i|) and (0.83±g;§f x 10" 3 , 
0.46^o;i5) for the power-law and log-normal models, respectively. The results for two models 
are consistent with each other. The power-law and log-normal models imply 1.9igi an d 
1.8l ;g times as many unbound or distant Jupiter-mass objects as the main sequence stars. 
a>3 is consistent with the values derived without planetary population indicating that brown 
dwarfs are 0.7 ± 0.3 times as common as main sequence stars. The numerical values of the 
models are summarized in supplementary Table 3. 
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Table 1. Microlensing Parameters of Short-Timescale Events. 
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Note. — N tE indicates the number of data points within to±^E and t Q , t-^, A mSu ^ and J s indicate the time of peak magnification, 
the Einstein radius crossing time, the maximum magnification, and the source star magnitude of the best fit models of the MOA 
data, respectively. JD'=JD-2450000. Uq and d m { n indicate the sourcelens impact parameter and minimum host star separation 
in units of the Einstein radii of the planetary mass lens, Re, and possible host star, R E *, respectively. The errors in £ E and u 
represent l-a limits. d min indicates the 2-a limits. MOA-ip-2, MOA-ip-3 and MOA-ip-10 were alerted as MOA-2007-BLG-144, 
MOA-2007-BLG-309 and MOA-2006-BLG-098 by the MOA real-time alert system dhttp://www.massey.ac.nz/~iabond/moa[ ). 
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This supplement to "Unbound or Distant Planetary Mass Population Detected by Grav- 
itational Microlensing" by Sumi et al. provides additional details of the analysis presented in 
this paper. The basics of gravitational microlensing is briefly introduced in Section [TJ The 
detailed data analysis and event selection methods are described in Section [21 In Section 
[H we discuss possible backgrounds that could conceivably result in false detections of short 
timescale microlensing event, and we show that the short timescale event sample does not 
have any significant contamination by false-positive detections. We discuss the simulations 
used for the detection efficiency calculations in Section HI In Section [51 we show that sys- 
tematic biases of the parameters produced by our analysis are small, and in Section [61 we 
present the details of the likelihood analysis that measures the substellar mass function. The 
likelihood analysis for other mass functions are also given in Section [7J Finally, in Section 
[S] we compare the frequency of unbound or distant planetary mass objects that of bound 
planets. 

1 Basic of the gravitational microlensing 

The gravitational microlensing method allows the detection of extremely faint or com- 
pletely dark objects when their gravitational fields act as a lens to magnify background 
source stars^. A gravitational lens is characterized by its Einstein ring radius, which is the 
radius of the ring image seen when the lens and source are perfectly aligned. In a microlens- 
ing event, the image separation is too small to resolve, and the observable feature is the 
changing total magnification due to the lens-source relative motion. The key parameter of a 
microlensing light curve is the Einstein radius crossing time given by 



where v t is the relative transverse velocity between the lens and source, Re is the Einstein 
radius, M is the lens mass, D\ and D s are the distance to the lens and the source from the 
observer, G is the gravitational constant and c is the speed of light. Although depends on 
the velocity and distance of the lens object, in addition to its mass, the distribution can be 
used to probe the mass function of the lens objects statistically. Thousands of microlensing 
events due to stars in the Galactic disk and bulge have been detected and used to study the 
stellar mass function and galactic structure^^£»^. 




, where 




(1) 
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2 Data Analysis and Event Selection 



The 2006-2007 MOA-II galactic bulge data consists of ~ 8250 images of each of the two 
most densely sampled fields (fields gb5 and gb9) and 1660-2980 images of each of the 20 other 
less densely sampled fields as shown in Table [TJ We observed these galactic bulge fields for 
about 8 months of each year from the end of February to the beginning of November. This 
observing strategy is designed to detect single lens events and short timescale anomalies in 
the light curves of stellar microlensing events due to planets orbiting the lens stars^Mi. The 
MOA images were reduced with MOA's implementation— of the difference image analysis 
(DIA) method— In the DIA method, a high quality, good seeing, reference image is 
subtracted from each observed image after transforming the reference image to give it the 
same seeing and photometric scaling as the observed image. This method generally provides 
more precise photometry in the very crowded Galactic bulge fields than straight PSF-fitting 
routines, such as DOPHOT 36 . This is, in part, due to the fact that the Galactic bulge fields 
are so crowded that virtually all the main sequence stars are not individually resolved. As 
described in Section HI the identification of a clear red clump giant population in the data is 
needed to match the observed MOA luminosity function to the much deeper Hubble Space 
Telescope luminosity function that describes the source stars. For the gb22 field and some 
fraction of other fields, totaling about 12% of the area, a clear red clump population could 
not be identified in the color magnitude diagram, and these regions were excluded from the 
analysis for this reason. This lack of a clear red clump giant feature in field gb22 is probably 
because it is relatively far from the center of the Galaxy. It is generally regions of very 
high interstellar extinction that prevented the identification of the red clump giant color 
magnitude diagram feature in portions of the other fields. 

The images were taken using the custom MOA-Red wide-band filter, which is equivalent 
with the sum of the standard Kron/Cousins R and /-bands. The instrumental magnitudes 
of the MOA reference images were calibrated to the Kron/Cousins /-band using OGLE-II 
photometry map of the Galactic bulged. The mean magnitude zero-point were estimated 
from the 30% of MOA-II fields which overlap with the OGLE-II map. We applied this mean 
zero-point to all fields. Here the error of the calibrated magnitude is estimated to be ~0.25 
mag from the standard deviation of zero-points in overlap fields. Although this calibration is 
approximate, this does not affect following analysis at all because the luminosity functions, 
which is the only part of our analysis requiring calibrated magnitudes, are calibrated by 
using red clump giants, as discussed in Section HI 

The problem of microlensing event selection has been addressed in previous microlensing 
analyses^^ 1 ^^^ 1 ^, anc [ we adopt event selection criteria that are similar to those used 
in a number of these previous papers. We employ three levels of selection criteria, or cuts, to 
distinguish genuine microlensing events from intrinsic variable stars and astrophysical and 
non-astrophysical artifacts. 

(1) Cut-0: On a subtracted image, variable objects can be seen as positive or negative 
point spread function (PSF) profiles depending on whether the target star is brighter 
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or fainter than in the reference image. We use a custom implementation of the IRAF 
task DAOFIND^i to detect these variable objects (like microlensing events), with the 
modification that both positive and negative PSF profiles are searched for simultane- 
ously. This algorithm finds difference image peaks with a signal to noise ratio (S/N) 
> 5 and then applies several additional criteria to avoid the detection of spurious vari- 
ations that are not associated with stellar variability, such as cosmic ray hits, satellite 
tracks and electrons leaked from the saturated images of bright stars. The positions of 
detected objects are checked against those obtained in previous reductions of the field. 
When no object is cross-referenced, the object is classified as new and added to the 
list of variable object positions. If an object has previously been detected, the number 
of detections for this object, N detect , is incremented. We require that N detect to be > 2. 
About 9.7 x 10 6 variable objects were detected at this stage of the analysis, including 
a number of image artifacts of various types. 

(2) Cut-1: Light curves of the candidates passing Cut-0 were then created using PSF fitting 
photometry on the difference images. Only objects with more than 500 photometric 
data points (from 1660-8250 images depending on the field) were used for further 
analysis. The photometric error bars were calibrated with constant stars in each sub- 
field. We place a 120-day moving window on each light curve and define the baseline 
flux to be the weighted average flux outside of that window. We require the baseline 
to have more than 10 data points and have Xout/dof < 3 for a constant flux fit. As 
other microlensing surveys have found, the error bar estimates from the photometry 
code provide only an approximate description of the photometric uncertainty for each 
measurement. The actual uncertainties depend upon the distribution of stars in the 
immediate vicinity of the target. We therefore define normalized uncertainties a\ = 

Uj-y/xout/dof, where Oi are the error bar of the measured flux Fi for the zth measurement. 
We then search for positive light curve "bumps" inside the 120-day window, with a 
"bump" defined as a brightening episode with more than 3 consecutive measurements 
with excess flux > 3cr' above the baseline flux, Fb ase - We define a statistic %3+ — 
Ej (Fi — -Fb aS e) /u'i summed over consecutive points with Fi — Fb asc > 3<r-. We require 
that a bump should have X3+ > 80 to pass this cut. 

(3) Cut-2: The candidate events that pass Cut-1 are fit with a 5-parameter microlensing 
model assuming a point source and a single lens object". The 5 parameters are the 
Einstein radius crossing time, t El the time of peak magnification, t , the source-lens 
impact parameter (in units of the Einstein radius), u , and the source and background 
fluxes, F s and F^. Candidate events are rejected when the Markov Chain Monte 
Carlo (MCMC) fitting method does not converge, because such events generally do 
not have well constrained parameters and may add a systematic bias in the observed 
t E distribution (see Section [5] for more details). We place constraints on the fit x 2 
for these microlensing fits. We require that x 2 /dof < 2 for both the full fit and for 
the regions near the peak with |t — to I < £e an d \t — to\ < 2tE- We require that to 
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be within the observation period, and we exclude very low-magnification events by 
requiring uq < 1 with an error bar of cr„ < 0.3. We also limit the best fit Einstein 
radius crossing times to the range 0.3 days < £e < 200 days with an error bar of 
< 50% and < 12 days, and we require that the source should have an (approximately 
calibrated) /-band magnitude brighter than 20.0. We require that the best fit source 
flux not exceed the cataloged flux on of an apparent source on the reference frame by 
more than a factor of 3. This allows for some "negative blending", which can occur if 
a source happens to be located at a local minima in the flux of unresolved stars. 

Of course, most of the cuts on the fit parameters exclude legitimate microlensing events, 
but we use them because they allow us to exclude the background of non-microlensing 
events with a relatively simple set of cuts. The cuts on are of particular interest 
because these will effect the shape of the t E distribution that is one of the major results 
of this analysis. The t E < 200 day cut (and the use of a 120-day window in Cut-0) 
excludes a large number of long events, which would be very important if we were 
trying to measure the high-mass end of the present-day mass function. This is not 
a problem for this analysis, which focuses on the low-mass end of the mass function. 
The removal of short events with £e < 0.3 days probably does not remove any events 
from this sample, but there are certainly real events of such short timescales. It will 
be important to remove this cut in future analyses of larger event samples that will 
measure the shape of the mass function of the unbound or distant planetary mass 
objects. 

Events with systematic residuals from the best fit model are also rejected, and this 
cut depends on the significance of the microlensing signal. We define the maximum 
number of consecutive measurements that are scattered from the best fit model with 
excess flux more than 2-a and 3-a by AT 2ct and N&j, respectively. We require that 
X 3+ > 70N 2a -500, and X 3+ > 45iV 3(T or N 3a < 2. 

All event selection criteria are summarized in Table [2] We have selected events which 
can be well fit with a point-source, single-lens model, and we have omitted binary lens 
events in this analysis by the cuts above, although some previously unknown planetary- 
mass binary lens events have been found in this analysis (Bennett et al. in preparation). 
Although we have identified more than a thousand microlensing candidates in this data set, 
only 474 high quality microlensing events have passed our relatively strict criteria. These 
strict criteria ensure that £g is well constrained for each event and that there is no significant 
contamination by mis-classified events, which might bias the £e distribution. This potential 
contamination by mis-classified events is discussed in the next Section. 

3 Background Events 

There are a variety of effects that can give signals that bear some resemblance to short 
timescale microlensing events, so it is important to ensure that our sample of events is 
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not contaminated by any of these effects. The potential background effects that we have 
considered are the following: 

(1) Cosmic-ray hits: It is rare for cosmic ray hits to give a signal with the same profile as the 
observed PSF. So, most cosmic rays are rejected due to their lack of a PSF-like shape. 
There was a previous claim of the detection of free-floating planetary mass objects 
by microlensing in the globular cluster M22 based on imaging by the Hubble Space 
Telescope (HST)—, but this claim was retracted within a few months 4 ™ as the putative 
microlensing events were artifacts caused by cosmic ray hits very close to the targets 
stars on two consecutive images. Our data set is not subject to such confusion for a 
number of reasons. First, unlike the MOA images, HST images are under-sampled, so 
that cosmic rays are much more difficult to reject based on their PSF shape in the HST 
images than in the MOA images. Second, the light curves presented in the M22 paper 
are maximally under-sampled, with only a pair of consecutive (cosmic-ray split) images 
that were brighter than the baseline brightness. In contrast, we explicitly require at 
least 4 consecutive observations that are > 3 a' above baseline and implicitly require 
much larger number with our cuts on the fit parameter error bars. Finally, the cosmic 
ray rate is much lower for ground-based images than for HST images. In fact, each of 
our 10 £g < 2 day events have at least 10 observations at significant magnification, so 
there is no chance of contamination by cosmic ray hits. 

(2) Fast Moving objects: Our pipeline conducts PSF photometry at the fixed position 
where the variable object was first detected on a difference image. If an object moves a 
significant fraction of PSF during several exposures, the photometry code will begin to 
underestimate its brightness due to the off-center PSF. With this photometry method, 
the linear movement of such an object can produce a light curve which has some 
resemblance to a microlensing light curve. Any small Solar System body, such as an 
asteroid or Kuiper belt object, with a proper motion of ~ lasec h , can mimic a short 
timescale microlensing event. Dust specs on the camera window can also produce 
similar light curves due to slight changes in the telescope pointing between different 
images of the same field. While light curves due to moving objects or dust specs can 
sometimes be reasonably well fit with a microlensing light curve model, they generally 
give unphysical light curve parameters, such as a large u Q value, a small t E value, 
and a fit source brightness, which to bright to be compatible with stars seen in the 
un-subtracted images. The moving objects in this data sets were all rejected by two 
independent methods. We have identified 3743 moving objects in this data set by 
direct measurements of the changing image positions on the difference images, and we 
found that all the light curves of these objects were rejected due to poor microlensing 
fits or unphysical microlensing model parameters. 

(3) Cataclysmic variables (CVs): The most important background for this analysis is a 
particular class of CVs that can have short timescale brightening episodes and repeat 
rarely enough so that they might have been missed despite the multi-year, high cadence 
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monitoring of the MOA-II bulge fields. Our primary means for identifying CVs is 
through the microlensing light curve fits. Most of the CVs can be identified from 
their highly asymmetric light curves, with a steep rise and a slow decline. However, 
for some CVs, this asymmetry is not seen due to gaps in light curve sampling or 
large photometric error bars. However, even poorly sampled CVs can generally be 
distinguished from microlensing light curves because a theoretical microlensing profile 
fit yields unphysical values such as a very large value of u and a sources with much 
brighter baseline fluxes than allowed by the reference images 

We identified 418 CVs in our sample from a visual inspection of light curves that failed 
the cuts on the quality of the microlensing fit or the upper limit on the brightness 
of the source star. The timescales of the brightening of these CVs range from hours 
to months. An important property of CVs is that they repeat^, and we have used 
this property to find an upper limit on the CV contamination of our sample. Using 
an expanded data set running through 2010, we found 208 CVs with more than one 
outburst. The average interval between microlensing events with the same source is 
about 100 thousand years if the lens objects are not physically related. The rate of 
repeating microlensing events due to wide binary lens or source systems is < 1 in a few 
hundred events, so very few, if any, of these events are likely to be actual microlensing 
events. Slight contamination of this CV sample with microlensing events is not a 
problem because these CV samples are only used for determining an upper limit of 
CVs which could pass our selection criteria. 

This upper limit is derived as follows. Because these CVs were identified by their 
multiple outbursts instead of the light curve shape, they can be used as a background 
sample to check how well we can reject CVs by their light-curve shapes. We decomposed 
the light curves of these multi-outburst CVs into 421 individual outbursts. The light 
curves of these outbursts were fit with microlensing models, and we found that none of 
them passed our selection criteria within 0.3 < £e < 2 days. With zero false positive in 
421 outbursts, Poisson statistics provides an upper limit of 2.3 at the 90% c.L. Then, 
the expected 90% c.l. upper limit of contamination among the 197 single-outburst CVs 
in our sample is 2.3x197/421=1.1, or 11% of our ten £e < 2 day events. 

After the MOA event selection was complete, we extended the light curves by adding 
the 2008-2010 MOA data and found that none of the ten t-^ < 2 day events have any 
other brightening episode. The MOA group also requested the OGLE light curves of 
these events. Seven of the ten short events were also observed by OGLE-III, the third 
phase of OGLE survey, and none of them have any other brightening in the 8-year 
OGLE-III light curves. Furthermore, the event MOA-ip-3 was also observed during 
OGLE-I and OGLE-II surveys and exhibits no further variations, which gives 16 years 
of coverage. Although the event MOA-ip-10 was not observed by OGLE-III, it was 
observed during OGLE-II survey in 1996-2000 and no brightening was seen. Finally, we 
checked the public database of the MACHO project, and found that the light curves of 
MOA-ip-3, MOA-ip-5, MOA-ip-6, and MOA-ip-8 had light curves with no photometric 
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variation observed during the 1993-1999 observing seasons. These data support our 
conclusion that these short events are not CVs. 

(4) Background supernovae: The dominant background in the microlensing survey toward 
the Large and Small Magellanic Clouds (LMC and SMC) was supernovae in galaxies 
at large distances behind the Magellanic Clouds, and methods to remove this back- 
ground from the microlensing sample by means of their asymmetric light curves were 
developed^^. 

The event rate of background supernovae in galactic bulge microlensing surveys is 
much smaller than that toward the Magellanic Clouds for the following reasons. The 
interstellar extinction toward the galactic bulge is much higher, Ay = 1 — 6 mag 47 , 
than towards the Magellanic Clouds with Ay ~ 0.3 mag—. Also, the LMC is ~ 6.25 
times more distant than the galactic bulge, and the SMC is even more distant, so the 
Magellanic Could source stars tend to be ~ 2 mag fainter. The ratio of supernovae to 
microlensing events is also dramatically reduced by the factor of ~ 100 ratio between 
the high bulge microlensing ratei 7 -* 1 ^ 1 ^ 1 ^ and the low rate towards the Magellanic 
Clouds^^. 

Although supernovae do not repeat like CVs, the shape of supernova light curves is 
somewhat similar to those of CVs. Both tend to be highly asymmetric with a steep 
rise and a slow decline. So that they can be rejected by the microlensing nts^ 4 ^^ in 
the same way as we have rejected CVs. 

Although supernovae do not pass our microlensing selection criteria, they could not 
contaminate our short timescale event sample even if they did, as they have timescales 
of ~30 days or longer. 

(5) Binary microlensing events: Binary lens systems generate small caustics when the 
separation of the lens masses is either much smaller or much larger than the Einstein 
radius. Wide binary events, however, can only mimic a short timescale single lens 
event when one of the lens masses has a low mass- 1 ^ 4 ^ So, wide binaries are not really 
a background, and we consider them in § |H] where we consider the distinction between 
very wide and free-floating planets. The 2-cx limits on the minimum separation to a 
host star are shown in Table 1 in the main paper. These limits are independent of the 
mass ratio and value for the host star in the limit where q <C 1. 

The situation with close binaries is different, because close binaries have two small 
caustics far from the center of mass that can give rise short duration lensing events even 
when the masses of the lens objects are large. We have fit all ten of the £e < 2 day events 
with close binary models and found that for most of the short events, a close-binary 
model is strongly excluded, being disfavored over the single- lens model by A\ 2 > 30. 
The only exceptions are events with relatively poor sampling. For events MOA-ip-2 and 
MOA-ip-6, the close binary models are disfavored by A% 2 = 17.2 and 9.8, respectively. 
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These imply formal probabilities of < 0.02 % and < 0.8 % for these events to be caused 
by close binaries. 

Event MOA-ip-5 is the only event for which the close binary model cannot be rejected. 
The best fit model for this event is actually a planetary wide binary model, which has 
a x 2 value lower than the single lens x 2 value by Ax 2 = 27.8. However, this level of y 2 
improvement is likely to be due to systematic photometry errors spread over a large 
number of data points, and we do not consider this a large enough x 2 improvement 
to indicate the detection of a host star. The best close binary model also gets some 
X 2 improvement from regions of the light curve that are distant from the peak shown 
in Fig. [H and it has a x 2 value that is better than the best fit single lens light curve 
by Ax 2 = 11.9 with 4 more fit parameters (as the magnification of the single lens 
model is too low for finite source effects to be important). Formally, this would be 
significant at the 98% confidence level, but only Ax 2 = 5.5 of this difference comes 
from the region within 10 days of the light curve peak. A Ax 2 = 5.5 difference is only 
significant at the ~ 70% confidence level, so we do not consider the close binary model 
to be be significantly favored for this event. The light curve coverage is insufficient to 
distinguish between the £e < 2 day single lens and close binary models. 

We can compare the case of MOA-ip-5 to the short events with better sampling in 
order to estimate the probability that MOA-ip-5 is due to a close binary lens system 
and is not a short timescale event. This analysis has revealed the 9 other events listed 
in Table 1 in the main paper, 2 wide separation planetary (or brown dwarf) binary 
events and one close binary (Bennett et al. in preparation). There is one other short 
event that is a poor fit to a single lens event, but a worse fit to a close binary. (It is 
probably a CV.) Thus, there are ~ 11 short events that are best fit by single lens or 
wide binary models, but only one that are best fit by a close binary. So, it is unlikely 
that MOA-ip-5 is a close binary event. The probability that it is a close binary can be 
quantified in an analysis like that by Bennett—, but it is simplest to just assume that 
MOA-ip-5 is a single lens event. 

We can also constrain the short event contamination due to close binaries by using 
the observed binary fraction and distribution. About 30% of main sequence stars 
have binary companions^. Only the small caustics due to close binaries with sepa- 
ration 0.05 < d/R-E < 0.3 in Einstein radius unit can resemble the short single-lens 
microlensing events. For 0.05 > d/R E , the caustics become too small to be detected. 
For d/R-E > 0.3, the effect of the host star is observable. This range corresponds 
— 1 < log(a/AU) < for the typical lens systems, i.e., K or M-dwarfs. Only 6% of 
binaries have separation within this range-^™ The event probability and detection 
efficiency for short events, which depend only on magnification timescale, are equiva- 
lent for close binary events and for the single lens events due to unbound or distant 
planetary mass objects. So, we can directly compare the number of close binary and 
single lens systems to estimate the contamination by close binaries. The fraction of 
short events due to close binaries versus single lens events is 0.3 x 0.06 x 2/1.8 = 0.02. 
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There is an additional factor of two for the close binaries because of the two small 
caustics generated a close binary lens system. This is consistent with our conclusion 
above that the contamination of our sample of short events by close binary lens systems 
is negligible. 

(6) Microlensing by high velocity stars and Galactic halo stellar remnants. Stars with very 
high velocity could also produce short timescale events. The mean timescale of the 
events due to stars in the bulge is t# = Re/vi = Re/{Di^) ~ 20 days, where \i = v t /D\ 
is the lens-source relative proper motion. So, it would require stars with at least an 
order of magnitude higher velocity than the typical stars to explain the events with 
£e < 2 days. 

For lens stars with high velocity at fixed distance D\, t& is proportional to pT . The 
typical proper motion of stars in our bulge fields is \i ~6 masyr -1 , so lens stars with 
jj, > 60 masyr -1 can have £e <2 days. However the number of such high proper motion 
stars with \i > 60 masyr -1 are only 6 x 10 -5 of stars in the OGLE-II proper motion 
catalog in the galactic bulge fields^. The event rate per lens star is proportional to the 
Einstein radius, Re, times a transverse velocity, v t , so the event rate due to the high 
proper motion stars with a typical mass of 0.3 M & is a factor of a/O.3M / a/O.OO1M x 
60/6 = 170 higher than that of a Jupiter mass lens with typical kinematics. If we make 
the (unrealistic) assumption that all the high proper motion stars are high velocity stars 
in the bulge, then we can obtain the expected rate of £e < 2 day events due to these 
high velocity stars. This is simply the fraction of high proper motion stars times the 
difference in event rate, or 6 x 10 -5 x 170 = 0.01 per star. This compares directly to 
the inferred number, 1.8, of Jupiter mass lens objects from the main analysis. So, the 
background of high velocity stellar lenses is more than two orders of magnitude less 
than the signal. Moreover, a more straight forward interpretation of these high proper 
motion stars in the OGLE-II proper motion catalog is that they are just just nearby 
stars instead of high velocity stars. Microlensing events due to nearby high proper 
motion stars are already included in our simulations of the £e distribution, so they do 
not produce the observed population of short £e events. Thus, the above upper limit 
of 0.01 per star is a very conservative limit. 

The conclusion that these high proper motion stars are nearby stars with regular 
kinematics instead of more distant high velocity stars are the following: 

(a) The proper motion distribution in the OGLE-II catalog of the bulge follows a slope 
of log(iV) = — 3 log(/i)+ const., without any distinct features, for fi > 10 masyr -1 . This 
matches the expectation for a standard kinematical model and a uniform distribution 
of stars in spaced. This implies that most of these high proper motion stars are not 
high velocity stars, but are instead nearby stars within about 300 pc which correspond 
the scale height of the disk. The proper motion of disk stars at a distance of < 300 pc 
with velocity of v ~ 30kms -1 is > 20 mas yr . 

(b) Galactic bulge lens stars contribute about 70% of all events observed towards the 
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bulge, and the typical bulge star velocity is about 120 km s -1 . A high velocity star, 
with a velocity an order of magnitude higher than this, would have a velocity higher 
than the escape velocity from the Galaxy. Galactic bulge radial velocity observations 
found that no stars with radial velocity of v\ > 400 km s -1 in a sample of 3200 stars 
and no significant deviations from a Gaussian distribution 55 . This confirms that there 
is no high velocity population in the bulge that could explain these short events. 

These arguments indicate that there is no significant population of high velocity stars 
to explain our detection of short events. 

A Galactic halo population of ancient stellar remnants, such as white dwarfs, neutron 
stars or black holes have sometimes been suggested as a explanation of the microlensing 
events seen towards the LMC by the MACHO Collaboration—. If so, their velocities 
should be v ~200kms _1 . The typical £g value due to the stellar remnants in the halo 
is > 20 days, so to get events with t E < 2 days, the velocities must be an order of 
magnitude larger, v ~2000 kms" 1 . So, like the stars considered above, a Galactic halo 
population of stellar remnants cannot explain these short events due to the Galactic 
escape velocity constraint. 

A final test of the possibility of high velocity lenses could be made by the detection of 
finite source effects for short events with relatively high magnification. Microlensing 
events due to massive lenses (with large Re) and high velocities should show this effect 
much less frequently than events due to planetary mass lenses (with small Re) and more 
typical velocities. In the current sample, events ip-1 and ip-10, have magnification high 
enough so that they might possibly show finite source effects. A finite source model for 
event ip-1, yields an improvement of A% 2 = 1.6 for a finite source model, which would 
imply a relative proper motion of fi ~ 6mas/yr, which is quite typical of bulge lensing 
events. For event ip-10, the best finite source model improves x 2 by only A% 2 = 0.1, 
but the implied lower limit on the relative proper motion is only /i > 4mas/yr, so most 
of the range of expected values is allowed. 

We can expand this test with two short events from the MOA alert system that occurred 
too late to be in our statistical sample. Event MOA-2009-BLG-450 has £e = 1.2 days 
and a peak magnification of 85. It shows no finite source effects, which implies a 
relative proper motion lower limits of fi > 6mas/yr. MOA-2010-BLG-418 was one of 
the shortest events known, with t E = 0.41 days and a peak magnification of 23. It does 
show a 3-0" signal for a finite source effect corresponding to a relative proper motion of 
fi = 7mas/yr. Thus, the event most sensitive to finite source effects does show them, 
and the other events have limits on fi which are consistent with expectations for a 
planetary mass population. 
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4 Detection Efficiency Simulations 



To compare observed t E distribution with the model, $ tE , we estimated the detection 
efficiency of our experiment using Monte Carlo simulations following previous work in the 
field~. Artificial microlensing events were added at random positions in the observed images, 
using PSFs derived from nearby stars in each field. The parameters of these artificial events 
were uniformly generated at random in the following ranges for the impact parameter, Uq, 
time of peak magnification, to, Einstein radius crossing time, t E , and source magnitude, /: 
< m < 1.5, 2453824 < t < 2454420 JD, 0.1 < t E < 250 days, and 14.25 < J < 21.15 
mag. (The to range is the range of observations in this data set.) The source magnitudes 
were weighted by the combined Luminosity function (LF) from MOA and the Hubble Space 
Telescope (HST)—. This uses the MOA LF for bright stars and HST for faint stars down 
to / = 24. This combined LF is calibrated to the extinction Galactic bulge distance for 
each field using the position of red clump giant stars in the color magnitude diagram, as red 
clump giant stars serve as a good standard candle. 

Once the images with artificial events were created, we processed them with the same 
analysis pipeline and selection criteria used in the analysis of the actual data. We evaluated 
our detection efficiency as a function of t E m each field by simulating 20 million artificial 
events. The results are shown in Fig. [2J While the detection efficiencies drop sharply for 
decreasing t E values, the detection efficiency is still significant at t E ~ 1 day, thanks to the 
high cadence of our survey. Fig. [3] shows the t E distribution with the data corrected for the 
detection efficiencies. This is similar to Fig. 2 of the main paper, except that in Fig. 2, 
it is the predicted model t E distributions that have been corrected for detection efficiencies. 
The flattening of the t E distribution at t E ~ 2 days indicates a rising mass function very low 
(planetary) masses, because the intrinsic microlensing probability scales with the lens mass 
M as ~ VM. 

5 Tests for Systematic Biases 

When light curve coverage is poor or the photometry errors are large, the light curve 
models can become degenerate. A model with a short timescale (small t E ) and a large ito 
value can look quite similar to a longer event with a smaller uq value and a fainter source 
flux, F s . The degeneracy tends to be worse for faint source stars blended with brighter 
stars, and it can lead to a bias in the inferred event parameters. This bias has been well 
studied—. To avoid such systematic biases in the observed t E distribution, we have rejected 
microlensing candidates with fit parameters that are not well constrained as described in 
Section [2j Note that the errors in t E listed in Table 1 in the main paper include the correlation 
with other parameters. In order to check for possible systematic parameter biases and 
significant contamination of our event sample, we have conducted several consistency checks 
with the observed and simulated events. 

The distribution of the impact parameter, uq which is purely geometric, can be used 
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as a test that the sample is dominated by real microlensing events. Fig. H] compares the 
observed and simulated uq distributions for events in different £e ranges. We performed the 
Kolmogorov-Smirnov test, which yields probabilities of 77%, 76%, 51%, 97% and 74% for 
the ranges £e =0.3-2.0, 2.0-13, 13-69, 89-200 days and the full sample, respectively. All 
the sub-samples are consistent with the simulated distributions, so there is no evidence of 
significant contamination by intrinsic variable stars in the Mo distributions. 

We also compared the input (true) t^in values input into our event simulations to the 
^E,out values output by the light curve analysis. The results are shown in Fig. [51 This test 
shows that the tE.out values are a reasonably accurate representation of the input £E,m values, 
but that there is a small systematic offset that makes the mean £E,out systematically larger 
than the input (true) t^ n values, independent of the value of £E,in- This implies that the 
bias in our measured £e is too small to influence our conclusions. Furthermore, the sign of 
the small bias that we do see implies that we might be slightly underestimating the number 
of short £e events. 



6 Likelihood Analysis of the Substellar and Stellar Mass 
Functions 

The Einstein radius crossing time is the only observable in the regular single lens 
microlensing event and is given by Equation ([1]). Although the physical parameters of the 
lens and source are degenerate, a model £e distribution, $(£e), can be calculated using a 
Monte Carlo simulation for an assumed mass function with a standard galactic mass density 
and velocity model—. The same method have been applied in the previous studies^ 1 ^ of 
the stellar mass function in the Galactic bulge. 

We used a Bayesian likelihood analysis to determine the mass function model parameters 
for both the power-law and log-normal mass functions. For the stellar and brown dwarf 
mass functions this was done both with and without the events with t-^ < 2 day events, 
and we found that the differences in the implied stellar and brown dwarf mass functions 
was negligible, as shown in Table [3j Fig. [6] shows the likelihood distribution for the mass 
function slope in the brown dwarf regime, «3, for the power-law model. This analysis finds 
«3 = — 0.48i[!;29, which is slightly higher than, but consistent with, previous estimates of 
«3 = 0-0.3^^^. Note that this a 3 values does not change when we consider additional 
mass function models in the following section. 

Fig. [7] is a contour plot of the likelihood distribution for the log-normal mass function 
parameters M c and cr c . The contours indicate the 68% and 95% confidence intervals and the 
cross indicates the maximum of the two-dimensional likelihood distribution. The log-normal 
parameter values determined by this analysis are a mean mass M c = 0.12 ± 0.03 and a 
log-normal width of er c = 0.76 jlo;f|. These are consistent with the values of M c = 0.079t '; l2 e 
and a c = O^t^oi determined in the previous work 22 . 

The power-law and log-normal models indicate that the number of brown dwarfs is 
0.73±o; 2 g and 0. 701^30 (respectively) times the number of main sequence stars in the mass 
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range 0.08 < M/M & < 1.0, which we denote by iV*. These values are higher than the value 
of 0.2 for the model (3) in Section mi The best model t E distributions, $(£e), are indicated 
as dashed lines in red and blue for the power-law and log-normal models, respectively, in 
Fig. 2 of the main paper and Fig. [3j As one can see, the both models represent data 
very well for t E > 2, but there is a significant excess at t E < 2 days. For t E < 2 days, 
we have detected 10 events compared to the predictions of 1.5 and 2.5 events without the 
planetary mass population for the power-law and log-normal models, respectively. The 
Poisson probability of detecting 10 events when only 1.5 or 2.5 events are predicted are 
4 x 10 -6 and 3 x 10~ 4 for the power-law and log-normal models, respectively. At i E < 1 days, 
we observe 3 events compared to predictions of 0.12 and 0.24 events for the planet-free 
power-law and log-normal models, respectively. The contributions of the best fit planetary 
mass distributions are indicated by the red and blue dotted curves for the power-law and 
log-normal mass function respectively, and the full distributions including lens objects of all 
masses are given by the red and blue solid curves. Both models provide an excellent fit to 
the observed t E distribution. 

7 Likelihood Analysis with the Other Mass Functions 

In this section we present the likelihood analysis for mass functions that differ from 
the models (1) and (2) shown in the main paper. Section 17.11 presents the results for an 
alternative stellar and brown dwarf mass function that is discontinuous near the brown 
dwarf-stellar boundary, and in Section I7T21 we consider a planetary power-law mass function. 
This is likely to be more realistic than the 5-function model, but its parameters are difficult 
to measure with the current data set. 

7.1 Discontinuous Substellar Mass Function 

In the main paper, we used two models, (1) the power- law mass function (MF) measured 
in the galactic bulge 21 with a continuous brown dwarf MF that was fit to the observed 
t E distribution, and (2) a log-normal MF fitted to the data, with best fit parameters are 
consistent with the values found with previous work—. 

Here, we present the likelihood analysis with model (3), a discontinuous MF— which 
uses the universal Galactic IMF, ip, based on the average of various stellar clusters^ 2 - with 
«! = 2.3 for 0.5M Q < M, a 2 = 1.3 for 0.075M© < M < O.5M , and a 3 = 0.3 for 0.01M Q < 
M < 0.075M Q . But this MF is discontinuous at the hydrogen burning limit Mhbl = 
O.O75M , with a discontinuity of a factor of -Rhbl = ^ 'a 3 (^hbl) / 'ipa 3 (A^hbl) =0.2-0.3. We 
assume that stars heavier than 1M have evolved to become stellar remnants, as we have 
done with the other mass functions. The steeper slope, oti = 2.3, implies fewer massive 
stellar remnants. In the planetary regime, we use a 5-function MF just as with models 
(1) and (2). Fig. [8] shows the result of likelihood analysis for -Rhbl = 0.3. The resulting 
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parameters are M/M Q = 1.9iJ;| x 10 -3 and $ = 0.50to;io as shown in Table[3j In this model, 
planetary mass objects with mass of 1.9lJ;g Jupiter mass are 1.3lo]I times as frequent as main 
sequence stars. With a slightly stronger mass function discontinuity, -Rhbl = 0.2, we find 
similar values: M/M & = 2.2±\* x 10~ 3 and $ = 0.5lloio- Fi g- E] presents the best fit t E 
distribution for the case of -Rhbl = 0.3. One can see that this model generates a narrower 
t E distribution. It has fewer long (>40 day) and short (<10 day) events compared to the 
peak (~20 days), than models (1) and (2) shown in the main paper. This is because this 
model has a steeper power law for high masses, O.5M < M, and fewer brown dwarfs. This 
model provides a worse fit to the data than models (1) and (2), with a maximum likelihood 
value that is a factor of 25 worse than model (1). However,the discontinuous MF model (3) 
does have one fewer model parameter that is fit to the data and it employs a mass function 
at the high end, which provided a worse fit to long timescale events in previous studies-^* 3 -. 
Since the relatively poor fit at the high mass end contributes to the poor likelihood value, 
we consider this model of the mass function of low-mass objects to be acceptable. The main 
conclusion, however, is that this mass function does not alter our conclusion that a new 
population of planetary mass objects is needed to explain the data, although it does imply 
a slightly higher mass for these objects. 

7.2 Power-law Planetary Mass Function 

In the main paper, we used a 5-function mass function for the planetary mass objects, but 
this is not likely to be realistic. Exoplanets are known to span a wide range of masses, and 
lower mass planets are much more easily ejected from their host stars in 2-body encounters. 
So, we have also considered a power-law form for the planetary part of the mass function, 
with an index of «pl for the mass range 1CT 5 < M/Mq < 0.01. This is assumed to be 
continuous with the power-law stellar and substellar mass function of MF model (1) in the 
main paper, at M — 0.01M Q . Our result does not depend on the lower mass limit 1O~ 5 M , 
because our sensitivity to events with M < 10 _4 M Q is quite small, due to the very low the 
detection efficiency at small t E . The power-law index, a^, for the brown dwarf mass regime, 
0.01 < M/Mq < 0.08, was fit with apL, simultaneously. 

Fig. [TUl shows the likelihood function for a 3 and «pl, and the resulting best fit param- 
eters are = 0. 49^^27 an d «pl = l-3loi- This value is consistent with our result for 
MF model (1) with and without 5-function planetary mass function (see Table |3]). The value 
apL, which is as steep as 0^2, is much steeper than a%, indicating that this planetary mass 
objects are separate population from the brown dwarfs. This model also predicts a larger, 
but poorly constrained, number of planetary mass objects per star with Npl/N* = 5.5 ^jf 3 , 
due to our poor sensitivity to very low-mass lenses. 

Fig. HH shows the best fit t E distribution compared with data and the best fit model 
(1) in the main paper for a comparison. Although this power- law planetary mass function 
model has a maximum likelihood that is a factor of 0.75 smaller than model (1), it has one 
fewer model parameter, so formally, it is a slightly better fit. This implies that the current 
data set is not able to constrain the shape of the mass function at sub- Jupiter masses. 
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8 Relative Frequencies of Unbound/Distant and Bound 
Planets 



Our sample includes three previously known events in which a planet and its host star 
are both detected: MOA-2007-BLG-192-^, OGLE-2007-BLG-368i, and OGLE-2007-BLG- 
349^. 

The short timescale events caused by planets bound to host stars with a separation 
of < 10 — 20 AU are detected as binary lens systems with a very small mass ratio. The 
events in our sample are detected as single lens events, because the planetary mass lenses 
are sufficiently isolated from any host star that the microlensing effect of the host star cannot 
be seen. These isolated planetary mass lenses can be either unbound or in distant orbits 
about a host star, and they can be distinguished from planetary mass lenses detected in 
binary events in two possible ways. The binary lens systems reveal the presence of a nearby 
host star by showing the effect of a planetary caustic in their light curves or by showing a 
low-magnification bump due to lensing by the host star— &§&L. We searched for the signature 
of host stars in our short events by fitting binary light curves to the data and found that 
none of the 10 selected events with t-& < 2 days shows any significant evidence of host stars 
(see bellow for their limits). However, our analysis has revealed three short events, which fail 
our selection criteria, that have clear binary lens caustic crossing features, as well as a very 
low amplitude signal due to lensing by the host star. Detailed modeling of these 3 events 
indicates wide binary models with lens mass ratios of q < 0.05 for two of these events and a 
close binary for the third (Bennett et al. in preparation). 

Thus, the total number of events in this sample in which both a planet and its host star 
are detected is 5. However, for two of these events, MOA-2007-BLG-192 and OGLE-2007- 
BLG-349, the planet is detected via light curve features due to the central caustic. This 
means that these planets have been detected only because of the gravitational effect of the 
host stars. For the other 3 planet plus host star events, the planet would have been detected 
without the influence of the host stars. Of course, there are likely to be some planets even 
closer to their host stars that are undetectable (by microlensing) because of the presence 
of the host star, so it is reasonable to presume that there are 5 planetary mass objects in 
binary lens systems in this sample that were detected with approximately the same efficiency 
as the planetary mass objects detected as isolated lens systems. This allows us to estimate 
the ratio of the new population of unbound or distant planetary mass objects to the planets 
found with host stars. The power-law mass function implies 11 isolated planetary mass lenses 
(that comprise our unbound or distant planetary mass sample) in our full event sample, so 
this ratio is 11/5 = 2.2. With the log-normal mass function, the number of events due to 
planetary mass objects is ~ 11, but 3 of these are from the low- mass tail of the brown dwarf 
mass function. This leaves 8 unbound or distant planetary mass objects for a ratio to the 
planets with detectable host stars of 8/5 = 1.6. 

The frequency of planets per star was measured by microlensing— for the planet-to-star 
mass-ratio interval —4.5 < log q < — 2 and separation a ~2-8 AU, to be d 2 N p \(d log qd log s) _1 = 



24 



0.36±0.15. This indicates N p \ ~0.9 for this q range and interval a = 1-10 AU. Radial velocity 
(RV) surveys have found that 11.3% of stars have planets with M smi > 10M e with periods 
less than 50 days^. The RV planets have periods that go down to about 3 days before the 
numbers really drop. So, we consider the RV survey to cover the period range 3 < P < 50 
day. If we assume a typical host star mass of 0.5 solar masses, then periods of P = 3 and 50 
days imply semi-major axes of a = 0.032 and 0.21 AU, respectively, corresponding a range 
of 0.81 in log(a). Assuming a uniform planet distribution in log(a), we have a 0.21 bound 
planets per star between a = 0.21 AU and 1 AU where the microlensing numbers take over, 
corresponding a range of 0.68 in log(a). If we assume a planet distribution in period of 
dA^pi(dlogP) -1 oc P - 26 as measured by RV for G dwarfs^ 8 , we have a 0.29 bound planets 
per star. Thus, we have a total of 1.13-1.19 bound planets per star in the semi-major axis 
range a = 0.03-10 AU. This implies 1.5-1.6 times as many planetary mass objects that serve 
as isolated lens systems to planets known to orbit stars. This result is nearly identical to 
the microlensing-only data based argument above. 

By necessity, we define our sample of isolated planetary mass lens objects as ones who 
don't have a host star that can be detected by microlensing, but this does not necessarily 
imply that host stars do not exist, which is why we interpret them as unbound or distant 
planets. The lower limit on the possible separation of a host star depends primarily on the 
magnification of the observed short timescale event. These minimum host star separations, 
G^min, are given in terms of the Einstein radius of the hypothetical host star, and the values 
for each of our £e < 2 day events are listed in the last column of Table 1 in the main paper. 
These are limits on the projected separation, and the typical semi-major axis corresponding 
to a host star with a separation of d min of one of our detected unbound or distant planetary 
mass object is ~ (3 AU)<i min , and therefore the typical lower limit on the separation of a host 
star is ~ 14 AU, with a range of ~ 7-45 AU depending on the event. 

Our data do not indicate what fraction of this population of unbound or distant plan- 
etary mass objects actually have host stars outside the region of detectability, but we can 
use the limits from direct detection searches to estimate this fraction. The Gemini Planet 
Imager- has set upper limits on the number of stars with Jupiter-mass planets at semi- 
major axes of 10-500 AU, with the tightest limits of < 30% of stars with Jupiter-mass 
planets with semi-major axes in the range 50-250 AU. Based on these results, we estimated 
that the fraction of stars with a 1 Jupiter-mass planet with semi-major axes in the range 
10 AU < a < 500 AU is less than 40% assuming a uniform distribution in log a. Since the 
new unbound or distant planetary mass population is comprised of ~ 1.8 times as many 
Jupiter-mass planets as stars, this comparison with the direct detection limits suggests that 
at least 75% of these unbound or distant planetary mass objects are not bound to any host 
star. 

Note that our likelihood analysis of the mass function in the planetary mass range, 
shown in Fig. [BJ givens a 1-cr lower bound on the 5-function mass of the planetary mass 
objects of approximately a Saturn mass. The direct imaging results don't rule out Saturn 
mass planets. So, if the majority of the lens objects in the newly discovery planetary mass 
population have a mass of order that of Saturn or less, then the direct detection limits do 
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not apply, and the majority of this new sample can be planets bound to host stars in wide 
orbits with semi-major axes of a ^ 7-45 AU depending on the event. 
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Figure 1. — Light curves of the 8 short microlensing events with ts < 2days that are not 
shown in the main paper. MOA data are in black and OGLE data are in red with error bars 
indicating s.e.m. For each event, the upper panel shows the full two-year light curve, and 
the middle and lower panels show a zoom-in around the peak and the residual from the best 
fit model, respectively. Although the light curve coverage of event MOA-ip-5 is somewhat 
sparse, it has 69 data points during the peak (to ± ^e), the second largest number in our 
short event sample (see Table 1 in the main paper) and our light curve model fits support 
the conclusion that this is a microlensing event. 
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Figure 1. — continued. 
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Figure 2. — The detection efficiencies of our experiment as a function of for the source 
stars down to I = 20.0 mag. Solid and dashed lines indicate the mean, minimum and 
maximum efficiencies of all fields. 
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Figure 3. — Observed and theoretical distributions of the event timescale £e- The black 
histogram represents the observed 474 microlensing events. This is the same data that is 
shown in Fig. 2 of the main paper, but here it is the data and not the models that are 
corrected for detection efficiencies. The error bars indicate s.e.m.. The best-fit models 
with the (1) Power-law and (2) Log-normal mass functions are indicated in red and blue, 
respectively. For each model, dashed lines indicate models for stellar, stellar remnant and 
brown dwarf populations and dotted lines represent the Planetary- mass population (PL). 
Solid lines are the sum of these populations. 
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Figure 4. — Cumulative distribution of the impact parameters, uq, for the observed sample of 
474 microlensing events (histograms) split into different £e bins. Dashed red curves indicate 
the expected distribution from the simulation in each £e- The probability that the observed 
samples are compatible with the model are estimated with the Kolmogorov-Smirnov test, 
which yields probabilities, Pks, of 77%, 76%, 51%, 97% and 74% for the i E ranges £e =0.3- 
2.0, 2.0-13, 13-69, 89-200 days and all sample, respectively. The detection efficiency is lower 
for events with large uq (low magnification). For the longest t E bin, we expect relatively 
large S/N even for large Uq. However, we have set a window of 120-days in cut-1, which the 
estimates the S/N by comparing to the deviation outside of the window. For the long i E , 
the tails of events extend outside of the window and decrease the S/N. This decreases the 
detection efficiency for events large uq and t E - 
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Figure 5. — The output timescale, tE,out, as a function of input iE,m of artificial events in our 
event simulations. The number density of events in each bin is shown as a logarithmic gray- 
scale, as indicated on the right. The artificial evens were generated uniformly in magnitude 
and log(iE)- Only events which passed the selection criteria are presented, and we have 
not weighted by the Luminosity function or true iE.in distribution. The yellow dashed line 
indicates £E,m = ^E,out, and the green solid line represents the mean of t^- m as a function of 
a mean of tE.out- We found £E,in is systematically ~ 5% smaller than tE.out regardless of £e- 
The blue solid lines indicate 90% interval in each tE.out bin. 
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Figure 6.— Likelihood distribution for the Power-law index in brown dwarf regime 
without the planetary mass population, including only the events with £e > 2 days. The 
vertical lines and colored regions indicate the median (a^ = 0.48) and 68% (0:3 = 0.10 and 
0.77) and 95% (03 = —0.44 and 1.01) confidence intervals. 



38 



CV2 



LO 



o 




0.01 0.1 1 

M c (M ) 

Figure 7. — Likelihood contours for the mean mass, M c , and variance, a c of log-normal mass 
functions without the planetary mass population, including only the events with > 2 days. 
Two contours indicate the 68% and 95% confidence intervals. "+" indicate the maximum 
likelihood point. The median and 68% confidence intervals are M c = O.I2I003 an d 0c = 
0.76^o;ig and these are consistent with previous work—. 
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Figure 8. — Likelihood contours for the planetary mass function parameters: the fraction of 
all objects that are in the planetary mass population, $pl, and their mass, Mpl, assuming 
a 5-function mass function. The two contours in each color indicate the 68% and 95% 
confidence levels. "+" indicate the maximum likelihood points. The red and blue contours 
represent the Power-law (1) and log- normal (2) mass functions also shown Fig. 3 of the 
main paper. The green lines indicate the discontinuous mass function model (3)— which 
uses the universal stellar mass function and discontinuous brown dwarf mass function with 
a power-law index of a 3 = 0.3 and a scale of -Rhbl = 0.3. The planetary mass functions are 
consistent each other despite the different brown dwarf mass functions. 
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t E (day) 

Figure 9. — Observed and theoretical distributions in the timescale ts as in Fig. 2 of the 
main paper. The black histogram represents the observed 474 microlensing events with 
the error bars indicating s.e.m.. The green lines indicates the best-fit for the discontinuous 
brown dwarf mass function— model (3) with -Rhbl = 0.3. The red line indicates the 
model (1) Power-law mass function. In each model, dashed lines indicate models for stellar, 
stellar remnant and brown dwarf populations. Dotted lines represent the Planetary-mass 
population (PL). Solid lines are the sum of these populations. Although the discontinuous 
brown dwarf model is a poorer fit to the data than model (1), the fit on the low-mass end 
is reasonable, and the mass function for the planetary mass objects is consistent with the 
results for models (1) and (2), which are presented in the main paper. 
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«3 

Figure 10. — Likelihood contours for the Power-law indices in the brown dwarf regime, 
and apL in the planetary-mass regime. These brown dwarf and planetary mas regimes span 
the mass ranges 0.01 < M/M Q < 0.08 and 10~ 5 < M/M < 0.01, respectively. The 
two contours indicate the 68% and 95% confidence levels, and "+" indicate the maximum 
likelihood point. These results do not depend on the lower mass limit of 10~ 5 M Q as the 
sensitivity to masses lower than 10~ 4 M Q is poor due to the small detection efficiencies at low 
t-£. The «3 distribution is consistent with the distribution for model (1) with and without 
the ^-function planetary mass function (see Tabled]). The value «pl is much steeper than 
a 3 , indicating that this planetary mass objects are separate population from brown dwarfs. 
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Figure 11. — Observed and theoretical distributions of the timescale £e- The black histogram 
represents the observed 474 microlensing events with the error bars indicating s.e.m.. The 
purple line indicates the best-fit model (4) with a power-law stellar and substellar mass 
function, as in model (1) , but with a continuous power-law mass function for in the planetary 
mass regime, 1(T 5 < M/M < 0.08 (Power-PL), instead of 5-function mass function of model 
(1). The red line indicates the power-law mass function model (1) with a 5-function planetary 
mass function. For each model, dashed lines indicate models for stellar, stellar remnant and 
brown dwarf populations, and dotted lines represent the Planetary-mass population. Solid 
lines are the sum of these populations. This power-law planetary mass function model (3) has 
a slightly smaller likelihood value than model (1), but because it has one fewer parameter, 
it is slightly favored, formally. 
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Table 1. MOA-II Galactic bulge fields 



Field R. A. (2000) Dec. (2000) N {ramc 



gbl 


17:47:31. 


.41 


-34:14:31 


.09 


2,253 


gb2 


17:54:01. 


,41 


-34:29:31 


.09 


2,386 


Kb3 

o 


17:54:01. 


.41 


-32:44:31 


.09 


2,067 


e:b4 

O 


17:54:01. 


.41 


-30:59:31 


.09 


2,985 


£rb5 

o 


17:54:01. 


.41 


-29:14:31 


.09 


8,229 


e;b6 

o 


17:54:01. 


.41 


-27:29:31 


.09 


1,779 


gb7 


18:00:01 


.41 


-32:44:31 


.09 


1,970 


gb8 


18:00:01 


,41 


-30:59:31 


.09 


2,139 


gb9 


18:00:01 


.41 


-29:14:31 


.09 


8,301 


gblO 


18:00:01 


,41 


-27:29:31 


.09 


1,992 


gbll 


18:06:01 


.41 


-32:44:31 


.09 


2,004 


gbl2 


18:06:01 


,41 


-30:59:31 


.09 


1,790 


gbl3 


18:06:01 


.41 


-29:14:31 


.09 


1,811 


gbl4 


18:06:01 


.41 


-27:29:31 


.09 


1,770 


gbl5 


18:06:01 


.41 


-25:44:31 


.09 


1,952 


gbl6 


18:12:01. 


.41 


-29:14:31 


.09 


1,756 


gbl7 


18:12:01. 


.41 


-27:29:31 


.09 


1,792 


gbl8 


18:12:01. 


.41 


-25:44:31 


.09 


1,799 


gbl9 


18:18:01. 


.41 


-25:29:31 


.09 


1,704 


gb20 


18:18:01. 


,41 


-23:44:31 


.09 


1,679 


gb21 


18:18:01. 


.41 


-21:59:31 


.09 


1,659 


gb22 


18:36:25. 


.41 


-23:53:31 


.09 


1,782 



Note. - - iVf ramc indicates the number of ob- 
served frames (exposures) during 2006-2007, 
which are used in this analysis. Fields gb5 
and gb9 were observed with a 10 minute ca- 
dence, while others were observed with a 50 
minute cadence. 
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Table 2. Event Selection Criteria 



level criteria 


comments 


CUtO iV de tect > 3 


Number of frames in which the object is detected. 



cutl iV data > 500 

iVout > 10 

xL/dof < 3 

-^bump — 1 

X3+ = gj (Fj ~ Fbasc) K > 80 

cut2 fitting converged 
X 2 /dof < 2 
X?/dof < 2 
Xs/dof < 2 
0.3 < t E < 200 days 
cr t Jt E < 0.5 
cr tB < 12 days 
3824 < t < 4420 JD' 
m < 1 
a uo < 0.3 
I s < 20.0 

(F s — F cat )/F CSLt < 3 
X3+ > 70N 2u -500 

X3+ > 45iV 3a OR iV 3a < 2 



Number of data points 

Number of data points outside of the 120-day window 
X 2 outside of the 120-day window 
Number of bumps in the window, where a bump 
has > 3 consecutive points > 3cr f above baseline 
Total significance of consecutive points with > 3a' 



Fits never converge if parameters are degenerate 
X 2 for all data 
X 2 for \t\ < t-E 
X 2 for \t\ < 2t E 

Einstein radius crossing timescale 
Error in t E 
Error in £g 

Peak should be within observational period 
The minimum impact parameter 
Error in uq 

Apparent /-band source magnitude 

Source flux should not greatly exceed catalog flux 

Exclude systematic residuals (depending on 

total significance) 
same as above 



Note. - JD'=JD-2450000. o\ = ^i^Xout/dof. N 2a and N 3a represent the maximum number 
of consecutive measurements which are scattered from the best fit model with an excess flux of 
more than 2-a and 3-<r, respectively. 
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Table 3. Mass Function 





Mass 


Function 


parameter 


Fraction 




(M Q ) 




(M and a are in M ) 


to 



1 40.0 < M 

8.00 < M < 40.0 
1.00 < M < 8.00 
0.70 < M < 1.00 
0.08 < M < 0.70 
0.01 < M < 0.08 
0.01 < M < 0.08 
M = M PL 

2 40.0 < M 

8.00 < M < 40.0 
1.00 < M < 8.00 
0.08 < M < 1.00 
0.01 < M < 0.08 
0.00 < M < 0.01 
M = M PL 



3 40.0 < M 


Gaussian 


Black hole (M r = 5, a r = 1) 




0.00060 


8.00 < M < 40.0 


Gaussian 


Neutron star (M r = 1.35, a T - 


= 0.04) 


0.0061 


1.00 < M < 8.00 


Gaussian 


White dwarf (M r = 0.6, a r = 


0.16) 


0.097 


0.50 < M < 1.00 


Power- law 


cti = 2.3 




1.0 


0.075 < M < 0.50 


Power- law 


0:2 = 1.3 






0.01 < M < 0.075 


Power-law 


a 3 = 0.3, -Rhbl = 0.3 




0.19 


M = M PL 


5-function 


Mp L = 1.9t^xlO- 3 ,$ PL = 


0.5018:15 


1-318:1 


4 0.08 < M 




same as model (1) 






0.01 < M < 0.08 


Power-law** 


a 3 = 0.491^ w/ PL 




0-7318:^ 


10~ 5 < M < 0.01 


Power-law** 


«pl = 1-3±8;1 w/ PL 




r r+18.1 
°-°-4.3 



Note. - - # is the model ID number: (1) power-law+<5-function, (2) log-normal+<S-function, (3) discontin- 
uous power-law+(5-function and (4) power-law+power-law. Gaussian: dN/dM — exp[(M — M r ) 2 /(2af )]. Func- 
tion types are: Power-law: dN/dlogM — M x ~ a . Log-normal: dN/dlog M = exp[(log M - log M c ) 2 /(2a 2 )]. 
Fraction is the number of objects relative to the number of main sequence stars with 0.08 < M/Mq < 1, iV* 
(0.075 < M/Mq < 1 for the model 3). The number of stellar remnants is estimated by extending the upper 
main sequence Power-law a = 2.0 (a — 2.3 for the model 3) through this higher mass regime. *: 0:3 or (Mpl, 
<I>pl) are fit to events with £e > 2 days without the planetary mass (5-function. **: and (Mpl, $pl) («pl 
for the model 3) are fit simultaneously for the full sample. * * *: M c — 0.12 and a c = 0.76 are held fixed when 
fitting Mpl and $pl- The fraction of planetary mass objects in model (4) is large because it extends down to 
M = 1O _5 M0 where have no sensitivity. 
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Gaussian 

Gaussian 

Gaussian 

Power-law 

Power-law 

Power-law* 

Power-law* 

5-function* 



Black hole (M r = 5, a r = 1) 
Neutron star (M r = 1.35, a T = 0.04) 
White dwarf (M r = 0.6, a T = 0.16) 
ai = 2.0 
a2 = 1.3 

a 3 = 0.48l8i? w/o PL 
a 3 = 0.50l°:io w/ PL 
M PL 



l.ilil x io- 3 ,$ 



PL 



= 0.4918:11 



0.0031 
0.021 
0.18 
1.0 



070+0.22 
■ lo -0.19 

0.7418:3° 

1.9151 



Gaussian 

Gaussian 

Gaussian 

Log-normal* 

Log-normal* 

Log-normal* 

5-function*** 



Black hole (M r = 5, a r = 1) 
Neutron star (M r = 1.35, a T = 0.04) 
White dwarf (M r = 0.6, a T = 0.16) 

M c 
M c 

M c — u.±z,_ 003 , 

M PT , = 0.8318:5? x 10~ 3 , $ PT , = 0.46: 



n 1 O+0.03 

1 2+°- 03 

al2+ 0.03 



(Tr 



(Tr 



= o.76l8:S 
= o.76l8:S 

- 76+ - 27 

10~ 3 ,$PL = 



-0.17 
-0.15 



0.0031 
0.021 
0.18 
1.0 

7n+ - 19 

u - ' u -0.30 

0.17I8: 24 



1., 



2+1.7 
-0.8 
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